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Abstract 
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We investigate the thermodynamical properties of nuclear matter in heavy ion col- 
lisions, in particular, with respect to questions of thermodynamical instability, phase 
transitions and fragmentation. For this we analyze results of relativistic transport cal- 

■ culations of Au + Au collisons at intermediate energies, separately for spectator and 
\ participant matter. On one hand, we determine local thermodynamical variables from 

■ the analysis of the local momentum distribution; on the other, we analyze fragment 
£NJ \ energy spectra in a blast model scenario. We find that the spectator represents an 

instable, equilibrized fragmenting source, while in the participant no such common 
source can be identified. Our results compare well with experimental determinations 

■ of temperatures and flow velocities. 

-xr : 

Introduction 

m ; 

Heavy ion collisions in the energy range up to a few GeV per nucleon have been investi- 
gated to obtain information on the nuclear matter equation-of-state (EOS) in a domain 

■ where hadronic degrees of freedom are thought to be the appropriate description. In such 
collisions several aspects of the EOS can be explored: high density and temperature in the 
compression phase and low density and possible signals of liquid-gas phase transitions in 
the expansion and the final fragmentation phase. Much work has been performed on the 
first question both experimentally, see e.g. [[[[], and theoretically, e.g. ||, by investigating 
in particular flow observables but also particle production [||. The investigation is compli- 
cated by the fact that at this stage equilibrium is not generally reached in medium energy 



collisions and therefore non-equilibrium effects have to taken into account in extracting 
the high density EOS [§. 

In this contribution we discuss the information that can be obtained in the expansion 
phase. Since the NN-interaction is qualitatively similar to inter-molecular forces with 
short range repulsion and longer range attraction, a phase diagram similar to a van-der- 
Waals system is to be expected. Thus there should exist first order phase transitions in 
bulk nuclear matter; however, the question is how such phase transitions occur in finite 
systems and what signals can be observed [[|. The expanding system will finally have 
to decay into fragments of different sizes (called multi-fragmentation) and information is 
sought from experimental measurements of fragment distributions, spectra, correlations, 
etc. In the interpretation of fragment spectra often the concept of an equilibrated freeze- 
out configuration is used, which is the point in the evolution of the system, where the 
fragments are formed and essentially cease to interact (except for Coulomb interactions). 
An important question is whether such a configuration exists, whether it is in thermal 
equilibrium, characterized by a unique temperature, and in chemical equilibrium, such 
that all fragment species are formed at the same instant. 

In this work the approach to answer these questions is to perform transport calcula- 
tions evolving the system microscopically from the initial stage through the formation of 



fragments. Such calculations are based on non-equilibrium formulations Q and no concept 
of equilibrium, freeze-out, etc. is neccessary. However, one may analyze the calculation 
with respect to questions, like the degree of equilibrium, whether a temperature can be 
defined, whether the system enters a region of instability, etc. One can also analyze the 
final fragmenting state using the same procedures as in experiment. By comparing the 
results of these two approaches we hope to understand better the questions involved in 
phase transitions in the finite system of heavy ion collisions. 

Transport calculations are usually done in a semiclassical approximation on the mean 
field level. However, fluctuations and correlations beyond the mean field level are impor- 
tant in the thermodynamically instable situation of a phase transition. In this work the 
approach is used |?j to gauge the numerical fluctuations arising from the use of the test 
particle method in such a way as to simulate the physical fluctuations. To go beyond 
this approach fluctuations have been included in the Boltzmann-Langevin approach M, 
but mainly schematic calulations have been done with this method. Recently together 
with the Catania group we have proposed a new and practical way to include fluctua- 
tions in a physical way in transport calculations Q. This procedure is not used here but 
will be employed in the future in somewhat less empirical investigation of fragmentation 
phenomena. 

Thermodynamical analysis 

We base our investigation on relativistic transport calculations of the Boltzmann-Nord- 
heim-Vlasov type (2|, 0]. The equation of motion for the one-body distribution function 
f(x,p*) reads 

><9* + (PIF^ + m*d£m*) flP*] f(x,p*) = I coll . (1) 

with the kinetic momentum p*^ = p^ + and the effective mass m* = m + E s , where 
T,^ and E s are the scalar and vector self energies, respectively. I co ii is a Boltzmann-type 
collision term respecting the Pauli-principle. For the self energies we have adopted the 
non- linear parametrization NL2 Pj. In ref. [0] we compared this parametrization to more 
realistic non-equilibrium self energies based on Dirac-Brueckner calculations. With respect 
to the thermodynamical variables discussed here, we found no essential differences between 
the two models, and thus we use the simpler NL2 here. The calculations are performed in 
the relativistic Landau- Vlasov approach which was described in detail in ref. |f[o| ]. It uses 
Gaussian test particles in coordinate and momentum space and thus allows to construct 
locally a smooth momentum distribution. 

The local thermodynamical properties are obtained from an analysis of the momentum 
dependence of f(x,p*) at a given space-time point x. The local 4-current and the local 
invariant density po = ^/T^W is directly calculated from f(x,p*), furthermore the energy 
density e and the pressure P are obtained from the energy momentum tensor. For the 
determination of a local temperature the local momentum distribution is subjected to a 
fit in terms of covariant hot Fermi-Dirac distributions of the form |11] 

n(x, p, T) = p-i — T (2) 

1 + exp —{ii*—p*u fJ ')/T 



with the temperature T and the effective chemical potential ^*{T) . The local streaming 
four-velocity is given as = j^/po- Then the temperature T is the only fit parameter 
to be directly determined from the phase space distribution. This temperature is a local 
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Figure 1: Phase space evolution of a semi-central Au + Au collision at 600 AMeV and 
b=4.5 fm. The left column shows the density distributions at certain times, the two right 
columns the momentum distribution at the center of mass (upper and middle row) and in 
the center of the spectator (lower row). The rightmost column shows the fit to determine 
the local temperature. 



thermodynamic temperature, which in the following we denote as T\ oc . Expression (Q) 
is appropriate for a system in local equilibrium. We model the anisotropic momentum 
distributions in the earlier stages of a heavy ion collision by counter-streaming or 'colliding' 
nuclear matter 11, 12, [lj|, i.e., by a superposition of two Fermi distributions n^ 12 ) = 
+ n( 2 ) — 5n( 12 ) , where e)n( 12 ) = V nW • n( 2 ) guarantees the Pauli principle and provides 
a smooth transition to one equilibrated system. In [11, H it has been demonstrated that 



this ansatz allows a reliable description of the participant and spectator matter at each 
stage of the reaction. 

An example of such an analysis is shown in Fig. 1 for a semi-central Au + Au collision. 
The local momentum distribution is given for characteristic space-time points in the fireball 
(participant) and in the spectator region. Also the fit is shown that determines the local 
temperature T/ oc . The participant region is not in equilibrium at the early times shown, 
however, the spectator represents a well equilibrated system, as does also the participant 
at later times. 

Experimentally much of the information about the thermodynamical behaviour in 
heavy ion collisions originates from the analysis of fragment observables. The physical 
and practical procedure how to properly decribe fragment production is still very much 
debated |7|, M. Here we use the simplest algorithm, namely a coalescence model, as we 
have decribed it in ref. fl. In brief, we apply phase space coalescence, i.e. nucleons form a 
fragment, if their positions and momenta (xi,pi) satisfy \xi — XA < R c and \pi — Pf \ < P c . 
R C ,P C are parameters which are fitted to reproduce the observed mass distributions and 
thus guarantee a good overall description of the fragment multiplicities. 

Fragment kinetic energy spectra have been analyzed experimentally in the Siemens- 
Rassmussen or blast model [[lj, |], |l5| . In this model the kinetic energies are interpreted in 
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Figure 2: Fragment kinetic energy spectra for a central Au + Au collisions at 600 AMeV 
for several fragment masses Af. These are shown for the coalescence model (full points) 
and with an additional binding energy criterion (open points). The lines represent the 
blast model fits to determine the slope temperature. 



terms of a thermalized freeze-out configuration, characterized by a common temperature 
and a radial flow, i.e. by an isotropically expanding source. In this model the kinetic 
energies are given by 
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where p and E are the center of mass momentum and the total energy of the particle 
with mass m, respectively, and where a = j/3p/T. The flow profile n{(3) is usually well 
parametrized as a Fermi-type function ffl]; however, the results are not very different when 
using a single flow velocity, i.e. n{(5) ~ 5((3 — /?/). One then has two parameters in the fit, 
namely /3/ and the temperature parameter in eq.(§), which we call T s i ope . It is, of course, 
not obvious that T s i ope represents a thermodynamical temperature, in fact, here we aim 
to find its significance. The expression (||) has been applied to kinetic energy spectra of 
all fragment masses simultaneously, yielding a global T s [ ope (global) , or to each fragment 
mass separately, giving T s i ope (Af). Examples of fits to the generated fragment spectra are 
shown in Fig. 2 for a central Au + Au collision. Spectra are shown with the coalescence 
model and with the additional condition of a positive binding energy. It is found that 
the temperature T s i ope extracted is almost independent of this (however, the number of 
fragments decreases due to sequential decay). 



Spectator matter 

The spectator is that part of the system which has not collided with the other nucleus, 
but which is nevertheless excited due to the shearing-off of part of the nucleus and due 
to absorption of participant particles. We expect from Fig. 1 that it represents a well 
equilibrated piece of nuclear matter at finite temperature. In Fig. 3 we show the evolution 
with time of the local temperature and the density for the spectator in semi-central Au+Au 
collisions at various incident energies. It is seen that after the time when the spectator is 
fully developed the densities and temperatures are rather independent of incident energy 
and remain rather constant for several tens of fm/c, making it an ideal system in order to 
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Figure 3: Local temperature (top) and density (bottom) evolution in the spectator in 
semi-central Au + Au reactions (b=4.5 fm) at different beam energies as a function of 
time. 



study the thermodynamical evolution of low-density, finite temperature nuclear matter. 
In Fig. 4 we show evolution of pressure and density with time as a parameter. We see 
that after about 45 fm/c the effective compressibility K ~ dP/dp\T becomes negative 
indicating that the system enters a region of spinodal instability and will subsequently 
break up into fragments. At this time we find densities of about p ~ 0.4 — 0.5po an d 
T ~ 5 — 6 MeV, which is in good agreement with findings of the ALADIN group based on 



isotope thermometers [16, 17]. The figure also shows isotherms of the bulk matter EOS at 
temperatures in the range of the spectator temperature. It is seen that the finite system 
follows the bulk material rather closely. 
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Figure 4: Density-pressure trajectory for the spectator matter in a semi-central Au on 
Au reaction at 600 A. MeV as in Fig. 3. The dotted curves are the nuclear matter equation 
of state for T = 5 and 9 MeV (lower and upper curve, respectively). 



Recently the ALADIN group has also determined kinetic energy spectra of spectator 
fragments |l7j and has extracted slope temperatures using eq.(^). Corresponding spectra, 
generated in the coalescence model, were shown in Fig. 2 together with the blast model fits. 
In Fig. 5 we show the slope temperatures separately for the different fragment masses and 
also the local temperature for comparison. For the nucleons we have T s i ope = (7.3 ± 3.5) 
MeV which is close to the local temperature T] oc = (5 — 6) MeV (see Fig. 3), as one would 
expect. On the other hand the slope temperatures of the fragments are considerably 
higher than those of the nucleons saturating for Af > 3 around T s i ope ~ 17 MeV. The 
experimental values from ALADIN 17] are also shown in Fig. 5. It can be seen that 
the slope temperatures from the theoretical calculations and those from the data agree 
extremely well. 
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Figure 5: Spectator slope temperatures for different fragment masses Af for the reaction 
as in Fig. 5 (diamonds). Also shown is the nucleon local temperature (circle) and the 
temperature obtained from a statistical model (gray band, see text). 

We see that T s [ ope for fragments differs from that for nucleons and also from 7} oc . Sim- 
ilar differences were seen by the ALADIN collaboration between T s i ope from the fragment 
spectra and temperatures determined from isotope ratios Ti so . The difference has been in- 
terpreted in ref . |l7]] in terms of the Goldhaber model [^, [l9| . When a system of fermions 
of given density and temperature breaks up the fragment momenta are approximately 
given by the sum of the momenta of the nucleons before the decay which leads to energy 
distributions which resemble Maxwellians and thus contribute to the slope temperature. 
To check this argument in this case we therefore initialized statistically a system of the 
mass and temperature of the spectator, and subjected it to the same fragmentation pro- 
cedure (coalescence) and to the same fit by eq. (||) as we did for the heavy ion collision. 
The slope temperatures obtained from this statistical model are given in Fig. 5 as a band, 
which corresponds to initializations covering the range of values in Fig. 3. It is seen that 
the model qualitatively explains the increase in the slope temperature relative to the local 
temperature and the increase with fragment mass relative to that for nucleons. A similar 



conclusion was drawn in ref [17| using the results from ref. [19|. This shows that T s i ope 
is not a thermodynamic temperature but simulates a higher temperature. This effect has 
been called "contribution of Fermi motion to the temperature". 



Participant matter 



The participant zone in a heavy ion collision constitutes another limiting, but still simple 
case for the investigation of the thermodynamical behaviour of nuclear matter. In contrast 
to the spectator zone one expects a compression-decompression cycle and thus richer 
phenomena with respect to fragmentation. The situation becomes particularly simple if 
we look at central collisions of symmetric systems |l4], [ij ^Ofl . We find that in Au + Au 
at 600A MeV a very well developed radial flow pattern appears after about 20 fm/c and 
the pressure becomes isotropic at about 35 fm/c indicating equilibriation. The number of 
collisions drops to small values at about 40 fm/c. This condition we shall call (nucleon) 
freeze-out. Thus equilibration and freeze-out occur rather simultaneously. We find a 
density at this stage of about normal nuclear density and a (local) temperature of about 
Tioc ~ 15 MeV in the mid-plane of the reaction. 

We now also apply the blast model of eq. (|3|) to fragment spectra generated in the 
coalescence model at the end of the collision at about 90 fm/c. The results for the slope 
temperature T s [ ope and the mean velocity (3f are shown in Fig. 6 for a common fit to 
all fragments with Aj > 2 for different incident energies. These are compared to the 
corresponding values extracted by the EOS [14] and FOPI Jl| collaborations. Within the 



uncertainities of the description by blast model fits, where the parameters T s i ope and (3f 
are not completely independent, there is qualitative agreement between calculation and 
experiment. 
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Figure 6: Slope temperatures (top) and radial flow velocities (bottom) from blast model 
fits to fragment (Af > 1) energy spectra for central collisions for different beam energies. 
They are compared to data from 14, HJ]. 



As was done for the spectator we also apply the blast model separately for the different 
fragment masses At. This is shown in Fig. 7 at 600 AMeV in the left column. We observe 
that slope temperatures increase and flow velocities fall with fragment mass in contrast to 
the behaviour for the spectator fragments in Fig. 5 where T s [ ope was about constant. A 
similar behaviour has been seen experimentally at 1 A. GeV in ref. |l4], [l]] and in calcula- 



tions in ref. [21|. Thus this behaviour cannot be interpreted as fragments originating from 
a common freeze-out configuration, i.e. from a fragmenting source. To arrive at an inter- 
pretation we have shown on the right column of Fig. 7 the local temperatures and flow 
velocities for different times before the nucleon freeze-out, i.e. for t' = tf reeze _ out — t, with 
t freeze-out ~ 35/m/c. For fragment masses Af > 1 the slope temperatures and velocities 
behave qualitatively very similar to the local temperatures and flow velocities at earlier 
times . This would suggest to interpret the results on the lhs of Fig. 7 as signifying that 
heavier fragments originate at times earlier than the nucleon freeze-out. This may not be 
unreasonable since in order to make a heavier fragment one needs higher densities which 
occur at earlier times and hence higher temperatures. However, this does not neccessarily 
imply that the fragments are really formed at this time, since fragments could hardly sur- 
vive such high temperatures. But it could mean that these fragments carry information 
about this stage of the collision. Similar conclusions were drawn in calculations in the 



QMD model [22]. In any case it means that in the participant region fragments are not 
formed in a common equlibrated freeze-out configuration, and that in such a situation 
slope temperatures have to be interpreted with great caution. 
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Figure 7: Slope temperatures (upper row) and flow velocities (lower row) for the same 
reaction as in fig. 4 at i£& e am = 0-6 AGeV. In the left column for blast model fits for 
different fragment masses and also for Af > 1 and for all fragments; in the right column 
the local values from the momentum distribution at times before the freeze-out. 



Summary 

To summarize fragmentation phenomena in heavy ion collisions are studied as a means 
to explore the phase diagram of hadronic matter. Here we have limited the discussion to 
the final stages of a collision and to the question of phase transitions and fragmentation, 
and to determine the thermodynamical properties of the fragmenting source. One way 
to do this experimentally is to investigate fragment kinetic energy spectra. In theoretical 
simulations the thermodynamical state can be obtained locally in space and time from 



the phase space distribution. In this work we have compared this with the information 
obtained from the generated fragment spectra. We apply this method to the spectator and 
participant regions of relativistic Au + ^4?x-collisions. We find that the spectator represents 
a well developed, equilibrated and instable fragmenting source. The difference in temper- 
ature determined from the local momentum space (or experimentally from the isotope 
ratios) and from the kinetic energy spectra can be attributed to the Fermi motion in the 
fragmenting source as discussed in a Goldhaber model. In the participant region the local 
temperature at the nucleon freeze-out and the slope temperature from fragment spectra 
are different from those of the spectator. The slope temperatures rise with fragment mass 
which might indicate that the fragments are not formed in a common, equilibrated source. 

We thank the ALADIN collaboration, in particular W. Trautmann and C. Schwarz, for 
helpful discussions. This work was supported in part by the German ministry of education 
and research BMBF under grant no. 06LM868I and grant no. 06TU887. 
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